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Résolution d’équations non linéaires: Introduction 


Soit f une fonction numérique d’une variable réelle. 
On cherche les racines simples de l’équation 

(1) f(x) = 0 


• Isoler les racines, c’est à dire trouver un intervalle [a, b] 
dans lequel a est l’unique racine réelle de (1). 

« Trouver cet intervalle : théorème des valeurs 
intermédiaires : 

a f(a) * f(b) <0 f admet un nombre impair de racines 
a Si f(a) * f(b) >0 f admet un nombre pair de racines 
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Introduction 


On supposera donc désormais avoir trouvé un intervalle [a, b] 
où f admet une unique racine simple et on supposera que f est 
définie, continue, et autant de fois continûment dérivable que 
nécessaire. 
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Introduction 


Les algorithmes classiques que nous allons étudier sont les 
suivants: 

O Méthode de la bissection 
O Méthode de Newton-Raphson 
O Méthode de la sécante 
O Méthode du point fixe. 
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Méthode de la bissection. 


On suppose que que f est continue dans [a, b] et que 
f{a).f(b) < 0 
o on pose c = 

a si f(c) = 0 Alors p-c 
a si f{a) * f(c) < 0 on remplace b par c 
a sinon on remplace a par c, 

o on continue cette opération jusqu’à ce qu’on trouve p avec 
la précision demandée. 
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Algorithme de bissection (ou de dichotomie.) 


Trouver une approximation de la solution de f(x) = 0 dans 
[a, b]. 

On construit une suite d’intervalles ([a n ,b n ])n contenant la 
racine, 

on teste que a n ou b n est le milieu de l’intervalle [a n ~ i , b n ~ i]- 

Entrées : a, b, e et N 0 

Sortie : la valeur approchée p : f(p) = 0 
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Algorithme de bissection (ou de dichotomie.) 


O Si f(à) = 0 imprimer la solution est a. Si f(b ) = 0 imprimer 
la solution est b, aller à 1 0 

O si f(b) * f(à) > 0, imprimer (pas de changement de signe). 

Aller à 10 

O poser N = 1 

O Tant que N < N 0 .faire les étapes 5 à 8 
O poser p = ^ 

O Si f(p) = 0 ou ^ < e, imprimer p. Aller à 10 
O poser N = N + 1 

O Si f(a) * f(p) > 0, alors poser a = p, sinon poser b = p 
Q Imprimer après N 0 itérations l’approximation obtenue est p 
et l’erreur maximale est ^ 

© Fin Jo*<v 
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Méthode de Newton-Raphson: 


Le principe consiste à construire une suite (x n ) n , telle que x n+ i 
soit l’intersection de la tangente à la courbe de f au point 
(x n , f{x n )) avec l’axe horizontal. 

On a: 

f A = (x 0 , f{x o)), B = (xi,0)e axe(Ox) 

\ Ae\ B € D : y = ax + b 

donc 

f f(x o) = ax 0 + b f 3= f'(x o) 
l 0 = ax : + b ^ \ x, = x 0 - 
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Algorithme de Newton-Raphson 


Entrées: une approximation initiale p 0 
e (la précision désirée) 

N 0 (le nombre maximum d’itérations) 

Sortie: valeur approchée de p ou un message d’échec 
O N= 1 

O Tant que N < N 0 , faire les étapes 3 à 6. 

O Poser p = p 0 - ^ 

O Sij p - p 0 |< e alors imprimer p, aller à l’étape 8. 

Q Poser N = N + 1 . 

O Poser p 0 = P- 

O Imprimer la méthode a échoué après N itérations. 
O Fin. 
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Méthode de la sécante 


La méthode de Newton-Raphson suppose le calcul de f'(p). 

On remplaçe dans la méthode de Newton f'(p n ) par 

f{Pn)-f(Pn- 1 ) 

Pn ~ Pn - 1 

L’équation de la sécante s’écrit : 

s(x) = f(p„) + (x * Po) f (P">- f (P°-'î 

Pn ~ Pn - 1 

Si s(p n+ 1 ) = 0 , on en déduit: 


Pn + 1 


Pn ~ f{Pn) 


Pn ~ Pn - 1 
f{Pn)~f(Pn- 1 ) 


□ ► ► < 1 ► < 1 ► 


'OQ>0 


S. EL HAJJI et T. GHEMIRES 


Résolution d’équations non linéaires 



Résolution d’équations non linéaires 

Résolution de systèmes d’équations non linéaires 


Algorithme de la sécante: 


Trouver une solution de f(x ) = 0 
Entrées: deux approximations initiales p 0 et Pi 
e (la précision désirée) 

N 0 (le nombre maximum d’itérations) 

Sortie:la valeur approchée de p ou un message d’échec 
O poser N = 1 , q 0 = f{p 0 ), Qi = f{pi ) 

O Tant que N < No + 1 , faire les étapes 3 à 6 
O poser p = pi-q 1 

O Si | p - Pî |< e alors imprimer p, aller à l’étape 8 
O Poser N = N + 1 

O Poser p 0 = p^ , q 0 = q^ , P^ = p,q i = f(p) 

O Imprimer la méthode a échoué après N 0 itérations 
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Méthode du point fixe 


Nous pouvons observer que la méthode de Newton peut 
s’interpréter comme 

f(x) 

Pn + 1 = g{Pn) OÙ g{x) = X- {yA)- 

Si la fonction g(x) est continue et 

si l’algorithme converge (c.à.d. p n -> p ), 

Alors, puique p n+1 = g(p n ), on a p : p = g{p) ; 
on dit que p est un point fixe de g. 
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Algorithme du point fixe 


Trouver une solution de g(x) = x 

Entrées: une approximation initiale po 
e(la précision désirée) 

N 0 le nombre maximale d’itérations 
Sortie: valeur approchée de p ou un message d’échec 
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Convergence et ordre de convergence. 
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Convergence et ordre de convergence. 


Théorème: 


Considérons le segment S = [p 0 - a,p 0 + a] c D; 
si F est contractante sur S et si | ^(pg) — Po I — (1-/0 a, 
alors l’itération p n+1 = F(p n ) de point initial p 0 : converge vers 
l’unique point fixe p G S de F. 


Théorème: 


Si F est différentiable au voisinage d’un point fixe p et si 
| F(p) |< 1 alors : 

31/ voisinage de p tels que po g V et p n + 1 = F(p n ) converge 
vers p. 
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Ordre de convergence. 


Définition 

Considérons une 
e n = p n ~ P ■ 

Sl { e*!, } COnV( 

linéairement vers 
ordre. 

Si on a { , 

l ( e n-ir 

d’ordre k (k le pli 

suite {p n } convergeant vers p et posons 

erge, on dit que la suite p n converge 
; p ou encore que la méthode est du premier 

| converge, alors la convergence est dite 
js grand possible.) 
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Ordre de convergence: Exemple. 


La méthode de Newton est une méthode de type point fixe avec 


F(x) = x 


f(x) 

nxy 


Si p est racine simple de f(x) = 0, alors f'(x*) / 0 et 
il existe un voisinage V de p tel que pour tout po e V, 
la suite ( p n ) n converge vers p et l’ordre de convergence est 2. 


Pour déterminer l’ordre de convergence, on utilise la formule de 
Taylor en p : 

F(x) = F(p) + F'(p)(x - x*) + F"(flx) (x ~ P)2 . 
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Résolution de systèmes d’équations non linéaires 


Nous considérons le problème suivant : 

F : R n ->• R n , trouver p e R n tel que F(p) = 0. 

Nous allons pour cela étendre au cas de la dimension n > 1 
certains des algorithmes proposés dans les sections 
précédentes. 

Pour k > 0, et D ç R n , 

C k (D ) = Ensemble des fonctions k fois continûment 
différentiables de R n dans R r restreintes à D. 
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Résolution de systèmes d’équations non linéaires 


Nous supposerons que FeC 1 (D). 

J F (x ) est la matrice jacobienne associée à F et évaluée au 
point x — [x u ...,x n ] T de R n . 

(Jf(x)),i = (£(x)), i,J= 1 n. 

Pour une norme vectorielle donnée ||.||, 

B(p- R) = {yeR n :\\y- P \\<R}. 

Boule ouverte de rayon R et de centre p 
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La méthode de Newton et ses variantes 


On peut étendre la méthode de Newton au cas vectoriel : 

x(°) g R n étant donné x, 

pour k = 0,1, jusqu’à convergence: 

résoudre J F (x^)ôx^ = -F(x^), 

poser x (/<+1) = xW + 6x^ k \ 

On doit donc résoudre un système linéaire de matrice J F (x^) 
à chaque itération k. 
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Exemple 


Exemple 

Considérons le système non linéaire 


e4+4 _i=o 
e x f-4 -i=o 


qui admet pour unique solution p = (0, 0) r . 
Dans ce cas, 


F(x) = [e x i 2+x 2 - 1 , e* 2 x 2 - 1], 
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Pour ||<5 (/f) || 2 < 1CT 10 ( Sx comme test d’arrêt ) 

Si *(°) = [0.1 , 0.1] r 

on obtient en 26 itérations le couple [0.1* 1CT 8 ,0.13* 10~ 8 ] 7 
ce qui démontre une convergence assez rapide. 

Si x(°) = [10, 1 0] r , 229 itérations sont nécessaires pour obtenir 
une solution comparable à la précédente, 

la méthode diverge si x(°) = [20, 20 ] T . 

Le comportement est cependant très sensible au choix de la 
donnée initiale. 
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Théorème 


Soit F : R n R n et F e C 1 sur un ouvert convexe D de R n qui 
contient p. 

Supposons que J F \p ) existe et qu’il existe des constantes 
R,C et L telles que 


j-' ( p> 


< C 


et 

IMfOO-JfOOII <L\\x-y\\ Vx, y g fi(p; R) 

Il existe alors r > 0 tel que, pour tout x(°) g fi(p; R) 
la suite est définie de façon unique et converge vers p avec 




< CL 


Ak) _ y(*) 
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Démonstration. 


On va montrer par récurrence sur k que : 
x^ +1 ) g B(p; r) , avec r = min(R, 1/(2 CL)). 
Pour x(°) G B(p; r) 
la matrice inverse J F \x(°î) existe. 

On a 


V(P)( J f( x(0) - J f(P)\\ < U F KH) 
<CLr<l 




et on 


déduit du Théorème que J F 1 (x^ 0 )) existe car 


J F \xW) 



Jf'(P) 


1 - 


||J f (x(0))-J f (p)|| 

< 2 

Jf» 

O 

C\J 

VI 
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Démonstration - suite 


Par conséquent, x (1) est bien défini et 

x^-p = x(°> - p - Jp 1 {x^)[F{xW) - F(p)\ 

En mettant en facteur Jp 1 (x (0) ) dans le membre de droite et en 
prenant les normes, on obtient 


r0)_ 


< 


J F \x^) F(p) - F{xW) - J F {xW) 


p-x(°) 


< 2 C- 


p-x(°) 


On a majoré le reste de la série de Taylor de F. 
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Démonstration - suite 


Comme de plus x (0) g B(p; r), on a \\x* - x (0) || < 1 /(2 CL), 
d’où 


r0)_ 


1 

< - 

“ 2 


X y '’ - P 

Ce qui assure que x (1) g B(p; r) 


-*(°) 


On montre de manière analogue que si on suppose la relation 
vraie pour un certain k, alors elle est encore vraie pour k + 1 . 
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Remarques 


Si *(°) est assez proche de la solution p et si la matrice 
jacobienne est inversible. 

La méthode de Newton converge de manière quadratique. 


Il faut noter que la résolution du système linéaire peut s’avérer 
excessivement coûteuse quand n devient grand. 

De plus, la matrice ^(x^^peut être mal conditionnée, ce qui 
rend difficile l’obtention d’une solution précise. 

Pour ces raisons, plusieurs versions modifiées de la méthode 
de Newton sont proposées. 
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Remarque 


Si on note rW = F(xW) le résidu à l’étape k, 

La méthode de Newton peut être récrite sous la forme 

{ld - J G {x {k) )) (* (/<+1) - x<*>) = -r (/() 

avec G(x) = x - F{x) 
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Méthodes de Newton modifiées. 


Ce résultat ne donne pas d’indication constructive sur la 
manière de calculer les incréments h. 

En diminuant les h, on peut diminuer l’erreur de troncature 
commise dans Méthodes de type sécante. 
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Application Contractante 


On dit qu’une application G : D c R n R n est contractante sur 
l’ensemble D 0 c D 

s’il existe une constante a < 1 telle que 
l|G(*)-G(y)|| < a \\x - y || pour tout x, y dans D 0 . 

Propriété (théorème de l’application contractante): Si 

G : D c R n -)• R n est une contraction sur un ensemble fermé 
Do c D telle que G(x) g D 0 pour tout x g D 0 , alors G admet un 
unique point fixe dans D 0 
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Méthodes de Point Fixe 


Théorème du point Fixe : On suppose que 6:DcR"->fi" 
possède un point fixe x à l’intérieur de D et que G est 
continûment différentiable dans un voisinage de p 
On note J a la jacobienne de G et on suppose que le rayon 
spectral p(J G (p)) < 1 . 

Alors il existe un voisinage S de x tel que S c D et, pour tout 
*(°) g S, la suite définie avant demeure dans D et converge 
vers x . 
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Exemple 


Considérons le système non linéaire 


F(x) = [xf + xf - 1 , 2x 1 + x 2 - I] 7 = (0, 0) 7 
dont les solutions sont x\ = (0, 1 ) T et xf = (4/5, — 3/5) r 


Utilisons pour le résoudre deux méthodes de point fixe 
respectivement définies par les fonctions d’itération 


1-x 2 


1~*2 


Gi(x) 



2 
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suite 


On peut vérifier que G-\ ( x *) = x, pour / = 1,2; 

les deux méthodes sont convergentes dans un voisinage de 
leur point fixe respectif car 

7).j G2 (x|)=(? 7) 

— fï 

~ V 3 


et donc p{ J G1 (x*)) = 0 et p(J G 2 (x 2 *)) 
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En exécutant le Programme avec une tolérance de 10~ 10 

Si x(°) = [-0.9, 0.9 ] T , le schéma converge vers x en 9 
itérations. 

Si *(°) = [0.9, 0.9 ] T , le schéma converge vers x en 115 
itérations 

Cette différence de comportement entre les deux suites 
s’explique par la différence entre les rayons spectraux des 
matrices d’itération correspondantes. • 
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Remarque: 


La méthode de Newton peut etre vue comme une méthode de 
point fixe associée à la fonction 

G n {x) = x-J f \x)F{x). 
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